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Q : Abstract 

^ ■ In this paper, we investigate a kind of special quintom model, which is made of a quintessence 

field (pi and a phantom field 02, and the potential function has the form of V{4>i — (j)^). This 
M^ \ kind of quintom fields can be separated into two kinds: the hessence model, which has the 

(•~^ ' state of (j)\> (j)2, and the hantom model with the state 0f < (j)^- We discuss the evolution of 

^O . these models in the w-w'plane {lo is the state equation of the dark energy, and uj' is its time 

derivative in unites of Hubble time), and find that according tow > —1 or < —1, and the 



^ 
^ 



Q ' potential of the quintom being climbed up or rolled down, the uj-uj' plane can be divided into 

^«0 . four parts. The late time attractor solution, if existing, is always quintessence-like or A-like 

^^ \ for hessence field, so the Big Rip doesn't exist. But for hantom field, its late time attractor 

1^ I solution can be phantom-like or A-like, and sometimes, the Big Rip is unavoidable. Then 

I ' we consider two special cases: one is the hessence field with an exponential potential, and 

2 ■ the other is with a power law potential. We investigate their evolution in the uj-uj' plane. 

'^ \ We also develop a theoretical method of constructing the hessence potential function directly 

^ • from the effective equation of state function uj{z). We apply our method to five kinds of 

parametrizations of equation of state parameter, where uj crossing —1 can exist, and find they 

k>>( I all can be realized. At last, we discuss the evolution of the perturbations of the quintom field, 

Vh ' and find the perturbations of the quintom 5q and the metric <1> are all finite even if at the 

state of a; = — 1 and uj' ^ 0. 
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1. Introduction 

Recent observations on the Type la Supernova (SNIa) [1], Cosmic Microwave Background Radiation 
(CMB) [2] and Large Scale Structure (LSS) [3] all suggest that the Universe mainly consists of dark energy 
(73%), dark matter (23%) and baryon matter (4%). How to understand the physics of the dark energy 
is an important issue, having the equation-of-state (EoS) uj < —1/3 and causing the recent accelerating 
expansion of the Universe. Several scenarios have been put forward as a possible explanation of it. A 
positive cosmological constant is the simplest candidate, but it needs the extreme fine tuning to account 
for the observed accelerating expansion of the Universe. As the alternative to the cosmological constant, 
a lot of dynamic models have been proposed, such as the quintessence models [4], which assume the dark 
energy is made of a light scalar field. These models can naturally get a state of —1 < cj < 0, but the 
state of a; < — 1 can't be realized, which makes many other possibilities have been considered such as the 
k-essence models [5] and the phantom models [6], which have the non-standard kinetic terms [5]. Besides 
these, some other models such as the generalized Chaplygin gas (GCG) models [7], the vector field models 
[8, 9] also be studied by a lot of authors. Although these models achieve some success, some problems also 
exist. 

One essential to understand the nature of the dark energy is to detect the value and evolution of its 
EoS. The observational data shows that the cosmological constant is a good candidate [10], which has the 
effective EoS oi to = —1. However, there are several evidences showing that the dark energy might evolve 
from a; > — 1 in the past to to < —1 today, and cross the critical state of w = — 1 in the intermediate 
redshift [11]. If such a result holds on with accumulation of observational data, this would be a great 
challenge to the current models of dark energy. It is obvious that the cosmological constant as a candidate 
will be excluded, and the dark energy must be dynamical. But the normal models such as the quintessence 
models, only give the state of — 1 < a; < 0. Although the k-essence models and the phantom models can 
get the state of a; < —1, the behavior of uj crossing —1 can not be realized [12]. So a lot of more complex 
models have been suggested to get around this [13]. Obviously, the most natural way is to consider a 
model with two real scalar fields. A lot of people have studied the so-called quintom model [14, 15], which 
is a hybrid of quintessence and phantom (so the name quintom) . Naively, We consider the action 

S = /AV^(-^ + £,. + £„), (1) 

where g is the determinant of the metric g^u, T^ is the Ricci scalar, Cde and Cm are the lagrangian densities 
of the dark energy and matter, respectively. The quintom dark energy has the lagrangian density 

Cde = Cq = \{d^<Plf - \{d^(t>2? - V^(01, ^2) , (2) 

where 0i and (/>2 are two real scalar fields and play the roles of quintessence and phantom field, respectively. 
Considering a spatially Flat-Robertson- Walker (FRW) Universe and assuming the scalar fields 0i and (/>2 
are homogeneous, one obtains the effective pressure and energy density of the quintom field 

PQ = \<Pi-\<Pl-y{<t^iA2) . (3) 



pQ = l<i>i-l4 + yi^i,<f>2) , (4) 

and the corresponding effective EoS is 

(jjQ = —. : . (5 j 

It is easily seen tliat ujq > — 1 when (p1 > 02 is satisfied, while ujq < — 1 when (p1 < </>2 is satisfied. It 
is obvious that the quintom is the simplest phenomenological model of the dark energy with ujq crossing 
— 1. The hybrid of (pi and (f)2 in the potential function makes the models being varied and complex, which 
prevents one from analyzing their general properties. So it is interesting to look for some kinds of quintom 
models with simple potentials. The cosmological evolution of the quintom model without direct coupling 
between (pi and 02 was studied in Ref. [15]. They showed that the transition from loq > — 1 to luq < — 1 or 
vice versa is possible in this type of models. But they also found that the late attractor solutions of these 
quintom fields are always phantom-like or A-like, which may lead to the Big Rip. The reason is simple: 
since the quintessence and phantom fields haven't direct coupling, the energy density of quintessence field 
(with the EoS uj > —1) decreases with time, but which is increasing for phantom field (with the EoS 
a; < — 1). So at last, the phantom field must be the dominant component, which may lead to the Big Rip. 
In this paper, we investigate another kind of quintom models with the potentials 

V{(Pi,^2) = V{(f>l-(t>'2). (6) 

In this kind of models, the fields (pi and (p2 couple by this potential function. Compared with the models 
in Ref. [15], these models are easily to discuss for their simple potentials. In Ref. [16], the authors found 
that this kind of models may be the local effective approximation of the D3-brane Universe. It is easily 
found that this model is equivalent with the dark energy made of a non-canonical complex scalar field 
^ = (pi + i(p2 in form with the lagrangian density 

1 



^de-^ 



id^<^y + {d^<p*r - v{<^' + <5*') , (7) 



which has been advised by Wei et al. in Ref. [17], where the authors found that this model can easily 
realize a state crossing the cosmological constant boundary. It is interesting that this model can avoid 
the difficulty of the Q-ball formation which gives trouble to the spintessence. Furthermore, by choosing a 
proper potential, this model can be described by a Chaplygin gas at late time. The authors also found that 
the Big Rip is avoided in the models with the exponential potential and the (inverse) power law potential 
in the special cases with (p\> (p2. 

The main task of this work is to investigate the general characters of this kind of quintom models 
with the potentials in Eq. (6). From the invariance under the transformation with hyperbolic function, we 
separate these models into two kinds: the hessence models with cp'i > (p2 and hantom models with (p1 < 02- 
By analyzing their evolution in the lu-u' plane, we find that if V is positive (negative), to' + 3{1 + uj){1 — uj) < 



(> 0) is satisfied. So the potential being climbed up or rolled down can be immediately judged by the 
value of the function u;' + 3(l +uj){l — uj). We also find that the hessence field always has a quintessence- 
like or A-like attractor solution, and the Big Rip is naturally avoided; but the hantom field always has 
phantom-like or A-like attractor solution, which may lead to the Big Rip. These characters can be seen 
clearly in two kinds of hessence models which we have investigated in this paper. After these, we study 
how to construct the potential of hessence directly from the effective EoS: uj{z). We apply our method 
to five kinds of parametrizations of the EoS parameter, where u crossing —1 can exist, and find they 
all can be easily realized in the hessence models. In the last part of this paper, we mainly discuss the 
evolution of perturbations of the quintom fields. By altering the forms of the evolutive equations, we find 
the divergence doesn't exist in these equations even if at the state of uj = —1 and uj' ^ 0. So the values of 
the perturbations are finite. 

The plan of this paper is as follows: in section 2, we review the evolutive equations of the quintom 
models, and separate them into two kinds: the hessence and the hantom models. In section 3, we investigate 
their evolution in the uj-uj' plane and analyze the general characters of their attractor solutions. In section 
4, we focus on two kinds of hessence models: one with the exponential potential and the other with the 
power law potential, and study their evolution in the uj-uj' plane. In section 5, we discuss the method to 
construct the potential of hessence field directly from the parameterized EoS and apply it to five kinds of 
parameterizations. In section 6, we investigate the perturbations of the quintom fields and their evolutive 
equations. At last, in section 7, we have a conclusion. 

We use the units h = c = 1 and adopt the metric convention as (-I-, —,—,—) throughout this paper. 

2. The Hessence and Hantom Models 

The quintom field here we consider has the lagrangian density as 

>Cq = 1(8,^1? - l{d,c^2f - V{<tA - <g) ■ (8) 

One can easily find that this lagrangian is invariant under the transformation 

4'! ~^ 4*1 cosh(ia) — 4>2 sinh(ia) , (9) 

</'2 — > — <^i sinh(ia) -|- (^2 cosh(ia) , (10) 

where a is constant. This property makes one can rewrite the lagrangian density (8) in another form 

Cq = ^he = \ [{^^.4? - <p\d^ef] - V{4) , (11) 

where we have introduced two new variables {(f), 6), i.e. 

01 = (^ cosh 6* , (/>2 = 0sinh0 , (12) 



which are defined by 



bl-(j)l, coth0 = (/)i/(/.2 . (13) 



These models are dubbed the "hessence" in Ref. [17]. But it is clear that this form requires an additional 
requirement, 0f > 02, on the quintom models. In another condition with (p1 < (p2, one can rewrite the 
lagrangian density in Eq. (8) in another form 

1 



^Q - Cha — 2 

here the variables (0, 6) are defined by 



{d,(py + ^\d^9Y -v{<p), (14) 



In this paper, we dub them "hantom". In the following discussion, we will find that the hessence and 
hantom have different properties, especially the late time attractor solutions. 

Hessence Models 

Let us restart our discussion from the action 

S = Jd''x./^(^-r^+Che + ^m^ , (16) 

where the lagrangian density of hessence field can be found in Eq. (11). Considering a spatially flat FRW 
Universe with metric 

ds^ = dt^ - a^{t)-/ijdx'dx^ , (17) 

where a{t) is the scale factor, and 7jj = (5*- denotes the flat background space. Assuming and 9 are 
homogeneous, from Eqs. (16) and (11), we obtain the equations of motion for (/) and 9 

4> + 3H(j) + 0^2 + dV/d(l) = , (18) 

(l)'^9 + (200 + 3H(I)'^)9 = , (19) 

where H = d/a is the Hubble parameter, an overdot denotes the derivatives with respect to cosmic time. 
The pressure and energy density of hessence field are 

Phe = ^ ('A' - 0'^') - ^(0) , Phe = l (0' - 0'^') + ^(0) , (20) 

respectively. Eq. (19) implies 

Q = a^(f)'^9 = const , (21) 

which is associated with the total conserved charge within the physical volume due to the internal symmetry 
[17], if consider the hessence as a non-canonical complex scalar field. It turns out 

A Q 



a^S^ 



(22) 



Substituting this into Eq. (18), we can rewrite the kinetic equation as 



which is equivalent to the energy conservation equation of the hessence ph(, + 3H(phe + Phe) = 0. The 
pressure, energy density and the EoS of the hessence are 

Phe = I'P^ - Sr72 - ^('^) ' Phe = \^^ " 1^ + ^(0) , (24) 



'J^he 



1 ■^' 



Q' 



v{<t>) 



Q' 



2a6(/.2 



+ V{<t>) 



(25) 



respectively. It is easily seen that u>he ^ —1 when 0^ > Q'^ / {a^ cfp') , while u>he ^ —1 when (fP' < Q'^ / {a^ (f)'^) . 
The transition occurs when 0^ = Q"^ / {a^ cfp') . In the case of Q = 0, the hessence becomes the quintessence 
model. If we define the effective potential 



Veff ^ V 



2a602 ' 



the kinetic equation (23) becomes 



+ 3i7(^ + dVeff/d(l) = . 



(26) 



(27) 



This is exactly the Klein-Gordon equation of quintessence field with the potential V{(j)) = 14 jj ((/>). So 
the field (j) will seek to roll towards in the minimum of its effective potential K;//, but which doesn't 
mean that </> will tend to roll towards in the minimum of its real potential V . This is the most important 
difference from the quintessence model. Then when the field rolls down to its potential, when it climbs up 
the potential, and how it influences the EoS of the hessence model? This is the main task of section 3. 
Hantom Models 



Now let's return to another case with (^f < (;^2- The action is 



S 



j d^Xyf^ 



IZ 



levrG 



I '~'ha ~r •'-m 



(28) 



where the lagrangian density of the hantom can be seen in Eq. (14), which follows the kinetic equations 

(29) 



+ 3i7^ + _gl_^ = 0, 



i6^3 



.3j,2/i 



where Q = a (f) 6, is the conserved charge. The pressure, energy density and EoS are 



1 



Pha 



+ 



2o6</>2 



vi<P), 



1 



Pha 



+ 



Q' 



2a^f 



+ vi<p) 



^ha 



+ 



Q' 
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V{4>) 



1 - 



+ 



Q' 



+ V{ct>) 



(30) 



respectively. It is easily seen that coha ^ —1 when (j)'^ < Q'^ / {a^ (j)"^) , and coha ^ —1 when 0^ > Q"^ / {a^ cfp') , 
which is inverse to the hessence models. In the case of Q = 0, the hantoni becomes the phantom field, 
which is also the origin of its name. If we define the effective potential 



the kinetic equation becomes 



"'"^"^S^p- <'i' 



• + 3H^ - dVeff/d4> = , (32) 



This is exactly the Klein-Gordon equation of phantom field with the potential V{(j)) = Veff{(p). So the 
field (f> will seek to climb up in the maximum of its effective potential V^ff, but which doesn't mean that 
(j) will tend to climb up in the maximum of its real potential V. Also, the discussion on the value of V, 
and its relation with to and lo' will be shown in the following section. 

There are two important characters we should notice: First, when encountering the condition (/)f = 02) 
one must return to the lagrangian with general form (8), which can't be discussed in hessence or hantom 
models with functions (p and 0, but here we don't discuss this condition in this paper. Second, in hantom, 
if (p is replaced by i(j), one finds the hantom becomes hessence model. So the hessence is enough to describe 
all kinds of quintom fields with potential form V{(pi — ^l); if the value of (j) being an imaginary number is 
allowed. 

3. The Evolution of EoS of the Quintom Fields 

Important observables to reveal the nature of dark energy are the EoS uj and its time derivative in units of 
Hubble time to' . The simplest model, cosmological constant, has the effective state of uj = —1 and uj' = 0, 
which corresponds to a fixed point in the uj-uj' plane. Generally, the dynamics model of dark energy shows 
a line in this plane, which describes the evolution of its EoS. Recently, it is shown that the simple scalar 
field models of dark energy occupy rather narrow regions in uj-uj' plane [18, 19]: the quintessence has the 
state of a; > —1, which only occupies the region of a;' > —3{1 — uj){l + uj) , and if the quintessence has the 
tracker behavior, the region decreases to be uj' > —{1 — uj){1+uj). A basic physics distinction in scalar field 
physics requires the precision on the dynamics to be of order cr{uj') ~ 2(1 + uj) < 0.1 [18]. The phantom 
field (uj < —1) occupies the region of uj' < —3(1 — uj){l + uj), and if the phantom has a tracker solution, 
the bound becomes uj' < 3uj{l — uj){l +uj). For the general k-essence model with the tracker behavior, 
the bound on u' is u' > ]-r^(l — i^)(l + uj). To confirm the quintom models, the crossing of cosmological 
constant must be found. So the dynamics of dark energy, especially at high redshift, is very important. 
The SNAP mission is expected to observe about 2000 SNIa each year, over a period of three years. Most of 
these SNIa are at the redshift z £ [0.2, 1.2]. The SNIa plus weak lensing methods conjoined can determine 
the present equation of state ratio, ujq, to 5%, and its time variation, uj', to 0.11 [20]. It has a powerful 
ability to differentiate the various dark energy models. 



In this section, we will extend the phase space analysis to the quintoni fields, where u) crossing —1 
exists. First we consider the hessence model, which has the kinetic equation 

O^ dV 



If w / — 1 is satisfied, one can define a function 

l+UJ 



l-LO 



w 


2a>i<j)'^ 




V 



(34) 



Then the kinetic equation (33) follows that 



I dlnx 1 V 

"^ 6 din a ~ 3HVI + UJ ' ^ ^ 

where a is the scale factor, and we have set the present scalar factor ag = 1. In the case of Q = 0, hessence 

becoming the normal quintessence field, the formula (35) leads to the relation [18] 



V /3ac2(1+l<j) / ldlnx\ 

where k^ = SvrG, and 17^ is the energy density of the quintessence field. The minus sign corresponds to 
(j) > {V < 0) and plus sign to the opposite. This can follow a constraint on to' , 

uj' > -3(1 - Lo) {I + uj) , (37) 

where and before to' = duj/dlna. This bound applies to a general class of quintessence field which 
monotonically rolls down the potential. 

Here we return to the general case of the formula (35) with Q ^ 0. Define a useful function c^ = p/p. 
If the matter is a kind of prefect liquid, this function is the adiabatic sound speed of this liquid. But for 
the scalar field, this function isn't a real speed. For the hessence field, it can be written as 

where p is the energy density of hessence, and to is its EoS. If the hessence returns to the quintessence 
field, one always has c^ < 1, which is for the field (/) monotonically rolls down (V < 0) the potential in the 
quintessence model with EoS uj > —1. And for phantom field with cj < — 1, c^ < 1 is also satisfied, since 
the phantom field climbs up the potential {V > 0). But here for the hessence model, c^ > 1 can exist, 
which only needs V{1 + uj) > 1. In these kinds of models, the function c^ isn't a physical speed. So the 
case of c^ > 1 is consistent with the Relativity Theory. Here this function reflects the evolutive direction 
of the potential of the hessence. In the case of cj > —1, i.e. the quintessence-like, c^ > 1 (c^ < 1) indicates 
that V > {V < 0), the field (p climbing up (rolling down) its potential, and in the case of oj < —1, i.e. 
the phantom-like, c^ > 1 (c^ < 1) indicates that V < {V > 0), the field (p rolling down (climbing up) its 
potential. Inserting this function into Eq. (35), one gets 

1 dlnx I — cf, ,„„, 

6 a In a 2V 



which can be rewritten as 



Using the relation of 



Eq. (40) fohows that 



1 (ilnx\ I ( _^ 2\ P 



' + 651^ /('-")=#>»■ («' 



d In X 2uj' 



dlna (l + a;)(l-w) 



(41) 



> T^-^T • (42) 



/// 

IV 



(1+..)(1-C.)(1-C2) (l-c2) 

So the uj-uj' plane is divided into four parts 

/: c2 < 1 & w > -1, w' > -3(l-tj)(l+a;) ; 

c\>\ k uj<-l, Lj' > -3{l-uj){l+uj) ; 

cl<l & w < -1, uj' < -3{1 - uj){l + uj) ; 

cl>l & uj>-l, uj' < -3(1 - u;){l + uj) . 
This can be seen clearly in Fig. [1]. From Eq. (38), one can easily find that y < is satisfied in Region / 
and //, the field rolling down the potential, and y > is satisfied in Region /// and IV, the field climbing 
up the potential. So from the value of the function uj' + 3(1 — uj){l + lo) being positive or negative, one 
can immediately judge how the field evolves at that time. This is one of the most important results in 
this section. In the case of Q = 0, the hessence returns to the quintessence field, and conditions of c^ < 1 
and UJ > —1 are always satisfied. So only Region / is allowed, and the bound of to' > —3(1 — uj){l + to) 
is satisfied, which is exactly same with Eq. (37). In the general case of Q 7^ 0, these four regions are all 
allowed. 

Now, let's focus on the issue: how the EoS crosses —1 in the uj-uj' plane? Assuming at some time, the 
hessence being in Region / with uj > —1 and c^ < 1, there are two ways for the field to run to the regions 
with UJ < —1: 

a) One is that, the field runs across the Critical Point {uj,uj') = (—1,0), and arrives in Region // or 
///. Unfortunately, this can't be realized in finite time. We use Taylor expansion of uj' at the state around 
the Critical Point, and keep the first two terms, 

duj' 



^'=^'\-^^ Ouj 



iuj + l)= b{uj + 1) , 
-1 



where 6 is a constant number. Using the definition of uj' , this equation yields that |u; + 1| = a!'. It is 
possible to get uj = —1 at a 7^ 0, only if 6 < is satisfied. In this condition, only if a ^ 00, uj = —1 can 
be got. So the hessence field can't cross the Critical Point in finite time. 

b) The other way is to cross the dot line and arrive in Region //. This is the only way for the hessence 
field to cross the state of w = —1. From Eq. (38), one can easily get 



2 w 



ci=uj- — . (43) 

3 1+c^ ^ ' 



When u! = —1 and to' ^ 0, c^ is divergent. In section 6, we will prove that this divergence doesn't yield 
the physical divergence of the perturbations of the dark energy. 

At last, we discuss the possible late time attractor solutions of the hessence field. There are three kinds 
of solutions: a) The hessence has not an oscillating EoS and the late time attractor is phantom- like with 
a; < — 1. Since lj' = is satisfied for the attractor, this solution must be in the Region ///; b) For the 
similar reason, if the attractor is quintessence-like with to > —1, it must stay in Region /; c) The other 
possibility is the A-like attractor. In this case, the hessence will run to the Critical Point {uj,uj') = (—1, 0). 
But the phantom-like attractor is difficult to realize. From the expressions of the pressure p and the 



energy density p of the hessence, one knows that only if ^ 



1^2 



< 



Q^ 



2a^4>- 



■, the attractor is phantom-like. If 



o 
> is satisfied for the attractor, the value of „^ ,^ will damp quickly with time, and at last, it is 



2aP$' 



always unavoidable to arrive at 



1 J,2 



> 



Q^ 



2^ ^ 2 >is^ ' '^hich is the quintessence-like result; on the other hand, if 
\^\ < is satisfied for the attractor, it is inevitable to run to the state of |(/)| = 0, which is forbidden by 
the definition of the hessence model. So in the hessence models, the Big Rip is avoided naturally, which 
has been discussed in some special examples in Ref. [17]. 

Now, let's discuss the hantom model in the similar way. The pressure and energy density of the 
hantom are 



p 



Q' 



'2*^ ^2a6(/,2 



V 



P 



'2*^ ^2a602 



+ V 



respectively. And the state equation is 



UJ 



1 -.^ 



Q' 



2a^(j)'' 



+ V 



Q' 
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V 



The kinetic equation is 



+ 3^^ + ;S^ 



dV 



(44) 



(45) 



(46) 



which follows that 



UJ 



> 



1 



(47) 



(l+c.)(l-^)(l-c2) 

where c^ is also defined by c^ = p/p, and the relation (38) is also satisfied. So in the hantom models, we 
also can divide the lo-lo' plane into the exactly same four parts as in Fig. [1]. One can easily find that 
y < is satisfied in Region / and //, the field rolling down the potential, and y > is satisfied in Region 
/// and IV, the field climbing up the potential. In the case of Q = 0, the hantom returns to the phantom 
field, and the conditions of c^ < 1 and a; < — 1 are always satisfied. So only Region /// is allowed, and the 
bound of uj' < —3(1 — uj){1 + uj) is satisfied. For the similar reason as before, the late time attractor of 
hantom can be phantom- like (Region ///) or A-like (Critical Point). The former can lead to the Big Rip 
at the late Universe. In the following sections, we only discuss the hessence models to avoid the Big Rip. 



4. Two kinds of Hessence Models 



In this section, we discuss two kinds of special potentials of the hessence fields. One is the model with an 
exponential potential 



where A is a dimensionless constant. The other is the model with a power law potential 



(48) 



(49) 



where n is a dimensionless positive constant. These two forms of potentials are the most popular models, 
which are discussed in the scalar dark energy. In this section, we will numerically solve the kinetic equation 
of the hessence with these two kinds of potential functions, and study the evolution of uj and uj' in detail 
to check the results we mentioned before. We focus on four special models: 

Model al: ^ > 0, V{(t)) = VQe-^""^ with A = 1.0, QVU^o) = 5, ^o = -1.4 ; 

4>o < 0, V{(P) = Voe-^'^'f' with A = 1.0, QV(Pt0o) = 0-5, ujq = -0.7 ; 

4>o > 0, y(,/.) = VoiK<l)r with n = 2, Qyipt<Pl) =5, c^o = -1-4 ; 



Model a2 
Model bl 
Model b2 



4>o < 0, Vi(l)) = yo('^<^)" with n = 2, QVU'^o) = 0.5, uq 



-0.7 



where (po is the field </> with the present value, pt is the present total energy density, and loq is the present 
EoS of the hessence. In all these models, we choose the present density parameters ^heo = 0.7 and 
f^mO = 0.3. The first two models have the exponential potentials, and the latter two ones have power law 
potentials. The present EoS in Models al and bl are phantom-like, and which are quintessence-like in 
Models a2 and b2. The EoS of the hessence is 



to 



1 V 



Q' 



2a602 



VicP) 



1 



Q' 



2a60^ 



+ vi<p) 



(50) 



In Fig. [2], we plot their evolution in uj-uj' plane, and find that, except the Model b2, the behavior of 
to crossing —1 exists in all these models. The Models al and a2 run to the same attractor solution of 
(uj,uj') = (—2/3,0), i.e. the quintessence-like solution, and the Models bl and b2 run to the same point of 
(uj,uj') = (—1,0), i.e. the A-like solution. These results are same with the conclusion in Ref. [17], where 
the authors found that the hessence models with exponential potentials have the stable attractor solutions 
with UJ = —1 + A^/3, and the models with power law potentials have the stable attractor solutions with 
u; = — 1. In these models, EoS crossing —1 obeys the second way: crossing the dot line (excluding the 
Critical Point [uj^lo') = (—1,0)). So the divergence of the function c^ exists. 

In the LO-uj' plane, the Models al and b2 stay in the region / and // at all times, so the condition 
V < holds for all time, and the fields roll down their potentials. But for the Models bl and a2, they 
run from the regions with to' < —3(1 -|- uj){l — uj) to the ones with uj' > —3(1 -|- cj)(l — uj), so the fields 
climb up at the beginning and then roll down the potentials. These can be seen clearly in Fig. [3], where 
we plot the evolution of function / = V / Hp with the scale factor. 



5. Construct the Potentials of the Hessence Fields 



Generally, the observed EoS of dark energy is a function of redshift z, and the function form of i^{z) 
depends on the parametrized model. Now, how to know the potential function from the observed uj{z)'I If 
realized, it will be a direct way to relate the observation and dark energy models. In Ref. [21], the authors 
suggested a theoretical method of constructing the quintessence potential T^(0) directly from the state 
function uj{z). Since lo < —1 can't be realized in the quintessence models, this method is effective only 
for the state of — 1 < lj < 1. But the recent observations mildly suggest that uj crossing —1 is existing. In 
this section, we will develop this method to construct the hessence potential V{(j)) directly from uj{z). We 
apply this method to five typical parametrizations. 

Consider the FRW Universe, which is dominated by the non-relativistic matter and a spatially homo- 
geneous hessence field (p. The Friedmann equation is 

^■2 



H = —{Pm + Phe) , 



(51) 



where pm and phe are the densities of matter and hessence, respectively. The pressure, energy density and 
EoS of the hessence field have been written in Eqs. (24) and (25), from which we have 



y{4>) = 2 '^-'^ ~ ^'^^^ p^"" ' 



(52) 



•^^ = Te;^ + (1 + '^he)phe ■ 



(53) 



These two equations relate the potential V and field (p to the only function ph^,. So the main task below 
is to solve the function form Phe{z) from the parametrized EoS u)hf,{z). The energy conservation equation 
of the hessence field is 



Phe +3H{phe +Phe) =0 



(54) 



which yields 



Phe{z) = P/ieOexp 



3 f\l+UJhe)dHl + z) 

Jo 



PheoEiz) , 



(55) 



where the subscript denotes the value of a quantity at the redshift z = (present). In term of uj^eiz), 
the potential can be written as a function of the redshift z: 



1 



V[<Piz)] = o(l - ^he)pheoE{z) . 



(56) 



With the help of pm = Pmo(l + z)'^ and Eq. (55), the Friedmann Eq. (51) becomes 



H{z) = Ho flrnoil + zf + ftheoEiz] 



1/2 



(57) 



where QmO and O,heo are the present relativity densities of matter and hessence, respectively. Using 
Eq. (53), we have 

1/2 

(58) 



d(j) 
dz 



-^ + {l+UJhe)phe 



{l + z)H{z) 

where the upper (lower) sign applies if </> > {(p < 0). Here we choose the lower sign to avoid the state of 
(j) = 0. It is helpful to define three dimensionless quantities (j), V and C 

C = K^Q'^/pheO ■ 



V ^ V/p, 



heO 



Eqs. (52) and (53) become 



Vs 



dz (1 + z) 



C{l+Zrcp-' + il+UJhe)E{z) 

ro(l + z)3 + E{z) 



1/2 



V[(l>] = ^{l-LOh,)E{z) , 



(59) 



(60) 



(61) 



where rg = Qmo/^heO is the energy density ratio of matter to hessence at present time. These two 
equations relate the hessence potential V{(p) to the EoS of the hessence ujfie{z). Given an effective ujfie{z), 
the construction Eqs. (60) and (61) allow us to construct the hessence potential V^(</'). Here we consider 
the construction process with the following five parametrization methods. 
The first model we consider is the EoS with constant value [22]: 

Model a: tOhe = ^o- If "^o > —1) a quintessence-like value, the construction of the potential can be 
easily realized with Q = 0, where the hessence returns to the quintessence field. This condition has been 
discussed in Ref. [21]. Here we consider another case with ujq = —1.2 < —1, a phantom-like value, and 
construct its potential function. In this case, the function E{z) has a simple form 



^(z) = (1 + z)3(i+"'«) . 

Then we consider three two-parameter models [23, 24, 25]: 
Model b: cohe = "^o + ^iz, and the function E{z) has the form 



E(z) = (l + z)3(i+^«-^i)e3"'i^ 



Model c: ujhe = <^o + "^i j^^j and the function E{z) has the form 

E{z) = (1 + z)3(l+«'o+«'i)e-3«'iTf^ ^ 
Model d: LO^f, = loq+ loi ln(l + z), and the function E{z) has the form 



(62) 



E{z) = {l + z) 



3{l+wo)+|wiln(l+2) 



(63) 



(64) 



(65) 



Inserting these into Eqs. (60) and (61), we can numerically evaluate the potential functions. In Fig. [4], 
we plot these parametrized EoS ujhe{z), where we have chosen the parameters uq = —1.2, loi = 0.5. We 



find that, in the Models b, c and d, uj crossing —1 exists. In Fig. [5], we plot the evolution of (j) with the 
redshift z, where we have chosen the parameters tq = 3/7, C = 5 and the present value (pQ = 1. One finds, 
with the increasing of redshift, the values of (j) monotonically increase in all these models. So the condition 
of (/) = is avoided. 

From Eq. (61), one can get the evolution of the potentials V{z) of the hessence fields, which have been 
shown in Fig. [6]. From this figure, one finds the potential is climbed up for all time in the Model a. 
And in the Models 6, c and d, the fields roll down the potentials at the higher redshift and climb up the 
potentials at the lower redshift. They all arrive at the lowest points of their potentials at the redshift, 
where w^g + 3(1 — a;/je)(l + oJhe) = is satisfied. These results are exactly what we expect: the Model a 
with ojhf. < — 1 and w^g = 0, stays in Region /// in the tOhe-^'he plane, so F > is satisfied for all time. 
But for the other models, we plot them in the ujhe-uj'f^^ plane in Fig. [7]. At the lower redshift, they all stay 
in Region ///, which makes y > is satisfied. And at the higher redshift, they all arrive in the Region /, 
where y < is satisfied. Combining Eqs. (60) and (61), the potential functions V{(p) can be got, which 
have been shown in Fig. [8]. One finds these potentials are not monotonic functions of </>, except the Model 
a, which are obviously different from the normal quintessence models [21]. 

Recently, a lot of authors have considered the dark energy with oscillating EoS [26]. They discussed 
that this kind of models give a naturally answer for the "coincidence problem" and "fine-tuning problem" of 
the dark energy. And in some models, it can naturally relate the early inflation and the recent accelerating 
expansion. The most interesting is that these models are likely to be marginally suggested by some 
observations [27]. The difficulty is that this kind of EoS is difficult to realize from the general potential 
function. Many periodic or nonmonotonic potentials have been put forward for quintessence fields, but 
rarely give rise to periodic ij{z). Here we consider a kind of oscillating parametrization: 
Model e: cuhe = wq + f^i sm{-p^). At the high redshift z ^ Zc, the oscillation of uJhe{z) disappears, and 
ujhe — "^o- The EoS is oscillating only when z < Zc- Here we choose parameters: loq = —0.7 loi = 0.5 
and Zc = 10, so Whe crossing —1 exists. And the present EoS is coheo = ~l-2, and w^^q = 0.5, which are 
same to the values in Models b, c and d. The observations mildly suggest that the EoS of dark energy 
crossed —1 at very recent, which had been regarded as the second cosmological coincidence problem. The 
parametrization in Model e gives a naturally answer for this problem. Using this LOfidz), we also can 
construct the potential of the hessence by applying the Eqs. (60) and (61), which has been plotted in 
Figs. [6] and [8]. We find although the potential V{(/)) shows an oscillating behavior, this oscillation is 
different from the simple sine or cosine function, which appears at the pseudo-Nambu-Goldstone boson 
(PNGB) field [28] with the potential V{(/)) = Vo[l + cos{(p/f)], where / is a (axion) symmetry energy 
scale. Here the potential V{(j)) of the hessence is an oscillating function with the increasing (or decreasing) 
amplitude. This suggests the method to build the potential of the scalar field dark energy, which can yield 
an oscillating EoS. 



6. The Perturbations of the Quintom Fields 

If the dark energy is a kind of dynamical field (or liquid) , it is necessary to consider the perturbations of 
it. These studies have been done for many kinds of dark energy models, such as the quintessence fields, 
the phantom fields, the k-essence fields, and so on. Some models have predicted too large perturbations 
of the dark energy or the background metric. For example, the GCG models can produce the oscillations 
or exponential blowup of the matter power spectrum, which is inconsistent with observations [29]; the 
Yang-Mills field models have the imaginary sound speed, which makes the perturbations of the intrinsic 
spatial curvature ^ increasing rapidly at recent epoch [9]. For many models, which allows the existence of 
EoS crossing —1, the perturbations of the dark energy may be divergent at the state of w^g = ~1 [30]. But 
whether or not, this divergence exists in our quintom models? In this section, we focus on this question 
by discussing the evolution of the perturbations of our quintom fields. In the conformal Newtonian gauge, 
the perturbed metric is given by 



ds^ = a^(T) 



1 + 2^)dT^ - (1 - 2^)dxUxi , (66) 



here we have used the conformal time r, which relates to the cosmic time by dt = adr. The gauge- 
invariant metric perturbation ^ is the Newtonian potential and <I> is the perturbation to the intrinsic 
spatial curvature. Always the background matters in the Universe are perfect fluids without anisotropic 
stress, which follows that <I> = ^. So there is only one perturbation function <I> in the metric (66). 

Using the notations of Ref. [31], the perturbations of the dark energy (including our quintom field) 
satisfy 

6'd, = -{!+ ujde){ede - 2>^') - m{cl - ^^)h, , (67) 

C = -n(\- 3u;de)ede + micl - OJde)0de + k^ ( T^^^de + 'A , (68) 



where 7i = a' /a, and the 'prime' denotes d/dr. Cg is the sound speed of the dark energy, which is defined by 
(?g = 5pde/5pde, and its value can't be larger than the light speed c. Here we have assumed zero anisotropic 
stress, which is the case for matter and simple dark energy models. The perturbation 5de is defined by 
^de = ^Pdel Pdei and Q^^ is the divergence of the velocity of the dark energy. We should point out that both 
scalar fields and fluids obey the same forms of these two equations, and the only difference comes from 
the term of c^. From Eq. (68), one can find that when c^ ^ cx), where Wf^e = ~1 and w^g 7^ are satisfied, 
one will get infinite Q'^^. Fortunately, this divergence can't lead to the divergence of (J^e- It is helpful to 
define another function 

7?de = (1 + ^de)Qde , (69) 

then Eqs. (67) and (68) become 

^'de = -^de + 3$'(1 + UJde) " 3^(0^ - u;de)Sde , (70) 



d'd^ = -n{l - 2,UJde)'&de + k^C%e + k\l + Ude)^ . (71) 

We find that the divergence at cjde = — 1 disappears. For the hessence field, we have 



6phe = 4>{d4>) + %^ + ^S4> - ^,p^ , (72) 



6p,. = m) + ^-'^S^-^4>'. (73) 

In the frame where the perturbations of the scalar field 6(/) and {6(/)) are negligible, the sound speed of the 
quintom becomes c^ ~ 1. So Eqs. (70) and (71) become 

She = -^he + 3$'(1 + UJhe) " 3^(1 - UJhe)She , (74) 

<, = -n{l - 3L0he)^he + k^She + k\l + LOhe)'^ . (75) 

In general the evolution of the perturbations can be numerically computed, which depends on the com- 
ponent in the Universe and the special quintom models. For a complete study on the perturbations, the 
evolution of the metric perturbation <I> should also been considered, which satisfies the equation [32, 30] 



$" + 3H ( 1 + ^) $' + ^A;2$ + 
PtJ Pt 



M) 



T-L^ + 2n' 
P't ' 



^ = AT^Ga'(6pt-^6pt] . (76) 



The pressure pt = X^iPii ^^"^ energy density pt = X^i Pi-, which should include the contributions of baryon, 
photon, neutrino, cold dark matter, and dark energy. Especially at late time of the Universe, the effect 
of dark energy is very important. Combining the Eqs. (74), (75) and (76), one can numerically solve the 
function ^{t) for special quintom models. Although we won't calculate them in this paper, it also can 
be found that the value of <I> is finite as if the total function p'^/ p't is finite. So even if c^ = p'del P'de ^^ 
divergent, $ is also finite if only p\l p\ isn't divergent. We should notice that, the Eqs. (74), (75) and (76) 
are also satisfied for the hantom models, where we only need to replace b^^ with b^a^ '^he with 'dha^ a-nd uj^e 
with i^ha- We remind that the perturbations <I> and ^ can directly influence the CMB anisotropy power 
spectrum by the integral- Sachs- Wolfe (ISW) effect [33], 

where Xlss is the conformal distance to the last scattering surface and ji the I'th spherical Bessel function. 
The ISW effect occurs because photons can gain energy as they travel through time-varying gravitational 
wells. One always solves the power spectrum Ci in the numerical methods [34], which can directly compare 
with the observations [2]. This is one of the most important way to study the dark energy models. 

7. Conclusion 



Xlss 




(77) 



Understanding the nature of dark energy is one of the most important issues in the modern cosmology. 
Until recently, the most effective way is to detect the EoS ujde and its time derivative oj'^^ by the observations 
on SNIa, CMB ,LSS and so on. There are mild evidences to show that tOde crossing —1 exists at the very 
low redshift, which makes the building of the dark energy models difficult. A simple quintessence, phantom 
or k-essence field is insufficient. Although the states of oj^e > — 1 and ujde < — 1 can be realized in these 
models, they all can't give a state of ujde crossing —1. A lot of more complex models have been suggested 
to account for this problem. Among them, the simplest one is the quintom model, which is a hybrid of 
quintessence and phantom fields. This kind of models have been discussed by a lot of authors. In this 
paper, we focused on a kind of special quintom, which has the potential T^((/>f — (p^), where (j)i and 02 are 
the quintessence and phantom fields, respectively. We investigated the general characters of this kind of 
models. 

The lagrangian densities of the quintom are invariant under the hypergeometric transformation between 
the fields (j)i and (/>2, which makes one can separate the quintom into two kinds: the hessence and the 
hantom. The former has the state of 0f > (/)2, and the latter has the state of 0f < (?i>2. We discussed their 
evolution in the ujq-uj'q plane, and found this plane can be divided into four parts according to the values 
of V and loq + 1 being larger or smaller than zero. The fact F > (< 0) denotes that the potential of the 
quintom is climbed up (rolled down), and the fact ujq > —1 (< —1) denotes the field is quintessence-like 
(phantom-like). From their kinetic equations, we found uj'q > —3(1 — iOQ){l + ^q) is satisfied for the case 
of y < 0, and uj'q < —3(1 — ujq){1 + ujq) is satisfied when y > 0, which directly relates the evolution 
of potential to the value of EoS ujq. We also found that, if the late time attractor solution exists, which 
is always quintessence-like or A-like for the hessence field, so the Big Rip is naturally avoided. But for 
hantom, this solution can be phantom-like or A-like. These characters are clearly shown in two hessence 
models with the exponential potential and power law potential. 

In this paper. We also developed a theoretical method of constructing the hessence potential directly 
from the observable EoS uj^eiz)- We applied our method to five kinds of parametrizations of EoS parameter, 
where iVhe crossing —1 can exist, and found they all can be realized in hessence models. Especially, the fifth 
model with the oscillating u;/je(z), we found although the potential V{(f)) shows an oscillating behavior, 
this oscillation is different from the simple sine or cosine function. Here the potential V{(j)) of the hessence 
is an oscillating function with the increasing (or decreasing) amplitude. In last part, we discussed the 
evolution of the perturbations of the quintom model, and found the perturbations of the quintom 6q and 
the metric <I> are all finite even if at the state of ujq = — 1 and uj'q ^ 0. We should notice that, in our 
discussion, we haven't considered the possible interaction between the quintom field and the background 
matter, which may show some new interesting characters [17]. 
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Figure 1: The uj-lij' plane of the hessence model. This plane is divided into four parts according to the 
values of c^ — 1 and a; + 1 being larger or smaller than zero. 
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Figure 2: The evolution of four hessence models in uj-uj' plane. The thin arrows denote the evolutive 
direction of uj and uj' with time. 
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Figure 3: The evolution of the function ^(0) in the four hessence models. Here the dimensionless function 
/(a) is defined by f = ^. 
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Figure 4: The EoS of five kinds of parametrization models. 
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Figure 5: Evolution of (j) with the redshift z of the hessence fields. 
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Figure 6: Evolution of potentials of the hessence models. 




Figure 7: The parametrization models in uj-uj' plane. The thin arrows denote the evolutive direction of uj 
and Lo' with the increasing redshift z. 





^.v^ 


Mod. a 


' 1 ' 1 ' 1 ' 1 


' 1 













Mod. b 


/■' ''•, 










Mod. c 


/ 










Mod.d 


/ 

\ 








2.0 


Mod. e 


/ 
/ 




— 










1.5 - 




/' ^'-^ 


,..'■■ 


_ 


o 

CD 




y 


-e- 






! 






> 


s 
1.0 - 


\^ 


1 

! 

1 ^ .^•'' 


^^' 


- 




^•^^^i:^--.,,^-'--.^-.-'' 






■ 


^ — ' 


^^^^^^ 




- 




0.5 - 

1 


1 





1 


- 



KCJ) 



Figure 8: Constructed hessence potentials F((/)). 



